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ABSTRACT 

Context. Very high dynamical range coronagraphs targeting direct exo-planet detection (10 9 ~ 10 10 contrast) at small angular separation (few 
A/D units) usually require an input wavefront quality on the order of ten thousandths of wavelength RMS. 

Aims. We propose a novel method based on a pre-optics setup that behaves partly as a low-efficiency coronagraph, and partly as a high- 
sensitivity wavefront aberration compensator (phase and amplitude). The combination of the two effects results in a highly accurate corrected 
wavefront. 

Methods. First, an (intensity-) unbalanced nulling interferometer (UNI) performs a rejection of part of the wavefront electric field. Then the 
recombined output wavefront has its input aberrations magnified. Because of the unbalanced recombination scheme, aberrations can be free of 
phase singular points (zeros) and can therefore be compensated by a downstream phase and amplitude correction (PAC) adaptive optics system, 
using two deformable mirrors. 

Results. In the image plane, the central star's peak intensity and the noise level of its speckled halo are reduced by the UNI-PAC combination: 
the output-corrected wavefront aberrations can be interpreted as an improved compensation of the initial (eventually already corrected) incident 
wavefront aberrations. 

Conclusions. The important conclusion is that not all the elements in the optical setup using UNI-PAC need to reach the Aj 10000 rms surface 
error quality. 
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1. Introduction 

Optical coronagraphy in space is one of the most use- 
ful methods to achieve high dynamic range observa- 
tions for the direct detection of extra-solar planets (e.g., 
ICoulter 20051 lAime & Vakili 2006b . A coronagraph can 
reduce a central star intensity and its diffracted halo 
light around the star where exo-planets would ap- 
pear. Several coronagraph designs using advanced focal 
plane masks (Guy on et al. 1999 IKuchner & Traub 2 002; 
IBaba & Murakami 20021 IRiaud e t al. 2003) have been re- 
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ported which reduce the light energy, sometimes to zero in 
theory, inside a re-imaged pupil plane called a Lyot stop, and 
also at a final image. Some techniques of nulling interferome- 
try, which has been mainly considered for mid-infrared long- 



baseline interferometers ( Mennesson, Leger & Olliver 2005 i 



has been applied to optical coronagraph in a single telescope 
using rotation-shearing interferometers ( Baudoz et al. 2 000; 
ITavrov et al. 200 5) or lateral-shearing interferometers with 
overlapped or separated sub-apertures dShao et al. 2 004 ; 
Nishikaw a et al. 2005b . Pupil function modification can re- 
duce only the halo intensity by using shaped or apodized 
pupils dKasdin et al. 20051 IGalicher et al. 20051 1. Wavefront 
phase control by an adaptive optics (AO) system can re- 
duce the halo intensity of a limited area, called "dark 
hole" (Malbet et al. 199Sb . Multiple-stage coronagraphs 
(lAime & Soummer 20041 ITolls et al. 20 05), and combinations 
of these methods, or pre-optics schemes have been proposed 
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to achieve a higher dynamic range (Nishika wa et al. 2 005 ; 
INishikawa & Murakami 20061 lAbe et al. 20061 ). 

In any case, the required dynamic range of 10 9 ~ 10 10 
for direct detection of earth-like planets by optical coro- 
nagraphs can be achieved only with a very high qual- 
ity wavefront of A/ 10000 rms and an intensity uniformity 
of 1/1000 rms JKuchner & Traub 20021 ILowman et al. 20041 >. 
Indeed, wavefront aberrations throughout the optical train 
produce a so-called speckle halo noise in the image plane 
(e.g. JBorde & Traub 20061) that prevents direct detection of the 
planet if these speckles are too bright. Manufacturing or pol- 
ishing accurate mirrors would be a good direction for small 
mirrors, while it would not be easy for a telescope primary 
to reach the required surface accuracy. Then an AO system 
would be used to achieve the required wavefront accuracy 
where both measurements and corrections of the wavefront are 
important. Giant or young planets brighter than 10 6 contrast 
are the target of ground-based telescopes in the infrared wave- 
length where wavefront correction by an AO system is also a 
key technology for use against atmospheric turbulence (e.g., 



UNI 



Mirror 1 



IMakidon et al~2 006; Serabyn et al. 2006|ITamura et al. 200"6t 



In an AO system, a deformable mirror (DM) is controlled 
by a signal from a wavefront sensor (WFS). Recently, control 
accuracy of a DM has been shown at a level of /1/10000 rms 
(e.g., IBrown et al. 20031 Evans et al. 20051 ). while some com- 
mercial wavefront sensors have been reported to achieve abso- 
lute wavefront measurement accuracy of /i/1000 rms by con- 
ventional Shack-Hartmann method (e.g., HASO HP 26 from 
Imagine Optic). The same goes for interferometric sensing, 
where repeatability on the order of A/l 0000 rms is commonly 
achieved (e.g., Zygo VeriFire AT). In the high dynamic range 
coronagraph regime, however, a non-common path error prob- 
lem, i.e., differences between optics of the WFS and the main 
path for the star and the planet images, should be considered in 
the pupil plane WFS such as the Shack-Hartmann sensor, the 



interferometer, or a curvature sensor (e.g., Guyon 2006). 

Solutions for the non-common path error problem are fo- 
cal plane WFSs flCodona & Angel 2004tlBorde & Traub 20061 
IGive'on et al. 20061 ) or speckle nulling, which features fo- 
cal plane speckle measurement and iterative wavefront con- 
trol. A dynamic range of 10 9 within the dark hole has been 
achieved after a few thousand iterations by an AO system 
(Balasubrama nian et al. 200"6b where the wavefront correction 
approaches the requirement. 

These focal-plane sensing approaches, however, meet other 
problems. When the speckle intensity level gets closer to the 
planet intensity, long exposure times are required to obtain 
better S/N measurements of the residual wavefront aberra- 
tions. In addition, the wavefront stability of the whole tele- 
scope system must be guaranteed during this time lapse. A 
multiple stage coronagraph reducing speckles with AOs at 
every stage dTolls et al. 20051) faces the same problem. In 
the case of multiple stages with approaches using nulling 
coronagraphs or nulling interferometers, no wavefront sens- 
ing has been considered before (e.g. INishikawa et al. 2 005 ; 
Nishikawa & Murakami 2006). But the same wavefront sens- 
ing problems would exist, i.e. the intensity in the pupil plane (at 
the Lyot stop) can be zero or very low, exhibiting a random in- 
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Fig. 1. Schematic of a possible optical layout for the UNI-PAC 
method (unbalanced nulling interferometer followed by a phase 
and amplitude correction adaptive optics). The various optical 
planes are those referred to in the text. 



tensity pattern where phase singularity at zero intensity points 
makes wavefront measurements and corrections difficult. Thus 
in the very high dynamic range optics for direct detection of 
exo-planets, maintaining wavefront quality remains an issue as 
well as the development of coronagraphs. 

We propose a pre-optics scheme for coronagraphs, using 
a combination of an unbalanced nulling interferometer (UNI) 
and a two-DM AO system (Fig[T|) for Phase and Amplitude 
Correction (PAC). Similarities can be found with previous 
studies, i.e. pre-optics by a balanced nulling interferometer 
( INishikawa et al. 20031 INishikawa & Murakami 20061 ). coron- 
agraphic pre-optics with non-zero amplitude distribution in the 
pupil after an interference and complex amplitude correction 
(Ab e et al. 20061 ). and a multiple-stage wavefront correction by 
commercial quality AO systems dTolls et al. 2005b . 

The pre-optics concept using a balanced nulling inter- 
ferometer in front of other coronagraphs was shown to 
be effective in absorbing part of the required dynamic 
range (INishikawa et al. 20031 INishikawa & Murakami 20061 ) 
but could not maintain sufficiently high wavefront aberrations 
after the interference to be corrected again by a secondary AO 
system. As a consequence, it puts some very drastic constraints 
on the precision and accuracy of the upstream AO system. 
Abe et al. (2006) showed that a pre-optics, that uses a pseudo- 
coronagraphic stage and a complex amplitude filter at the Lyot 
plane, can produce a flat and almost uniform wavefront, in 
spite of central obscuration and spider-arm patterns (i.e. large 
amplitude aberrations at the entrance pupil). We extended this 
concept in order to facilitate the phase and amplitude measure- 
ment/correction after a nulling stage, which is usually the bot- 
tleneck of high contrast coronagraphs. In the proposed concept, 
the UNI aims at simultaneously performing the absorption of 
part of the dynamic range, and the magnification of wavefront 
aberrations without phase singularities for later compensation 
with an AO system (i.e. the PAC), so that the star's speckled 
halo is effectively attenuated. 

This paper introduces the principle of the UNI-PAC con- 
cept, a novel method that avoids the non-common path error 
problem and the low intensity or phase singularity problem. 
The overview, formalization, and simulation of the concept are 
addressed in Sect.[2j Sect.[3] and Sect.HJ respectively, followed 
by a discussion in Sect.|5J Conclusions are drawn in Sect.|6] 
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2. Overview of the concept 

A rough illustration of the UNI-PAC process is as follows. 
Figure[T]shows a possible optical concept of the UNI-PAC sys- 
tem, and Fig.|2]depicts the associated effects on the amplitude 
and phase in the pupil plane at different points in the set-up. The 
UNI uses two wavefronts extracted from a collimated beam of 
a telescope. This beam can be compensated by an AO system, 
for instance, at a reference level of A I 1000 rms that we will use 
throughout this paper. 

In our optical scheme, the UNI stage is intended to produce 
two effects. Firstly, it reduces the intensity of the star image 
by, e.g., an order of 1/100, but has no effect on the speckle 
halo level. This corresponds to an amplitude reduction of the 
unaberrated electric field by a partial nulling in the pupil plane. 

A second effect of the UNI is the wavefront aberration mag- 
nification (for instance, from A/ 1 000 to A/ 1 00), which is in fact 
a consequence of the partial nulling, and one of the most impor- 
tant phenomenon in the UNI-PAC concept. It is consistent with 
that the difference between the attenuated star intensity and the 
speckle level becomes smaller. 

The next stage is the PAC, that is, a second AO system 
whose goal is to compensate for the magnified phase and am- 
plitude aberrations (from Aj 100 to A/ 1000 in our example). 
The corresponding effect in the image plane is the reduction 
of the speckled halo level. The PAC stage can effectively con- 
trol the wavefront because the UNI stage has maintained an 
almost uniform intensity level in the pupil plane (i.e. with no 
phase singularities). 

In our example, before the UNI stage, the initial in- 
tensity ratios of the star's peak intensity, the speckle halo, 
and the planet are 1 : 10 7 : 10~ 9 under a wavefront aber- 
ration of about A 1 1000 rms. After the UNI, the ratios be- 
come 10~ 2 : 10~ 7 : 10~ 9 . Here the starlight is reduced but the 
speckle is still there. Then, after the PAC stage they become 
10~ 2 : 10~ 9 : 10~ 9 . Thus the UNI-PAC has reduced both the 
star and its speckled halo intensity levels by a factor of 100. 

A downstream coronagraph can then be used to further de- 
crease the star intensity to detect the planet. Note that the coro- 
nagraph does not require Aj 10000 optics quality any more, but 
only about Aj 1000 which is enough to perform a 10~ 7 starlight 
attenuation, and keep the speckle level down to the planet level. 

We also recall that the aberration level was /1/1000 both 
at the UNI-PAC and the first AO system locations (which can 
also include the non-common path errors of the AO systems). 
Despite an overall optical quality of Aj 1000, we have reached 
a dynamic range of 10 9 , meaning that the wavefront quality is 
"virtually" improved beyond the AO system capabilities. A big 
advantage of using the UNI-PAC is this specification relaxation 
for the entire optics of the coronagraph, which is probably key 
to achieving the planet detection system in a cost-effective way. 

3. Formalism 

3.1. Definitions and preparations 

We consider the complex amplitude of an electric field for each 
point on a pupil plane whose coordinate is (x, y) defined by the 
sum of an unaberrated constant amplitude E° and a complex 
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Fig. 2. Schematic of wavefront phase and amplitude magnifi- 
cation in the process of the UNI-PAC method, where only the 
real part is indicated for the amplitude. In planes A, B, B', B" 
solid and dashed lines show wavefront 1 and 2, respectively. In 
planes UNI and PAC the combined wavefront is shown. 

aberration term s gathering both amplitude and phase aberra- 
tions, 



E(x,y) = E° + s(x,y). 



(1) 



We prefer to use this expression rather than the somewhat clas- 
sical E = E° ( 1 + a) e" 9 notation (where a and 9 are the ampli- 
tude and the phase aberrations, respectively). Indeed, the ra- 
tio s(x,y)/E°, "the normalized complex aberration", is a use- 
ful indicator of the wavefront quality in our algebraic deriva- 
tions. Also, under small wavefront aberration conditions, its 
real and imaginary parts are almost equal to a(x, y) and 0(x, y), 
respectively (the correspondence between these notations can 
be found in Appendix lAl. That way, we can treat "the normal- 
ized complex aberration" as the wavefront aberration. We can 
assume E° is a real number and under a very small wavefront 
aberration condition 



E° = (E(x,y)) 
(s(x,y)) = 0, 



(2) 



where ( ) is a spatial ensemble average within the pupil plane. 

We use notations for the intensity of each complex ampli- 
tude of electric field as 

7° 

2 ( 3 ) 



= |£°l 2 
I s (x,y) = \e(x,y)\ 



where | | is the absolute value for real numbers or the mod- 
ulus for complex numbers. I s can be observed as a residual 
intensity in the pupil plane when the unaberrated component is 
sufficiently attenuated. The average of the residual intensity is 
proportional to the average speckle background intensity in the 
image plane. 

In the following equations, the notations described above 
will include a suffix, indicating the location where the wave- 
front is evaluated, that is, one of six planes (A, B, B', B", UNI 
and then PAC planes) in the optical layout of the concept of 
Fig. Q] In Fig. [3] a vector representation of the complex ampli- 
tude is shown for each plane in the UNI-PAC. 
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Fig. 3. Complex amplitude at a point (x, y) in the pupil plane 
under the processes of the UNI-PAC method. In each panel, an 
arrow from the origin indicates the unaberrated amplitude, E°, 
and the complex aberration component, e, is shown by other 
small arrow(s), where the distribution area of the complex am- 
plitude is drawn by a dashed-lined ellipse. Traditional phase 
aberration, 6, is indicated by an arc. The upper four panels in- 
dicate the complex amplitude of beam A, beam B with ND 
filter before n phase shift, UNI, and finally PAC planes respec- 
tively from top to bottom, in E A units. The two lower panels 
are similar to the last two frames of the upper panel, except 
that the vectors have been normalized by the factor gE A , where 
g-Q.2. The normalized aberration is magnified by 1 jg by the 
UNI (UNI output plane), and compensated by a factor h(x,y) 
at the PAC output plane. 



We consider two wavefronts E A (x,y) and E B (x,y) at the 
input of our system generated from a collimated beam of a sin- 
gle telescope. These wavefronts are corrected in advance by 
an AO system before the UNI stage, but aberrations (including 



non-common path errors) still remain written, 

( E A (x,y) = E°. +s A (x,y) 

\ (4) 
[E B (x,y) = E°+s B (x,y). 

Any AO system has a limited number of actuators, hence usu- 
ally limiting the correction to lower spatial frequency ranges. 
This is a common situation of AO systems, but fortunately the 
low spatial frequency region corresponds to the relevant inner 
area in the image plane where the exo-planets would be ob- 
served. In this paper, we assume that a rather ordinary qual- 
ity of /l/1000rms level for the AO system is used on the tele- 
scope side before our concept, and another AO of the same 
quality is again adopted in the latter half of our concept (the 
PAC stage). We also consider how the combination of the un- 
balanced nulling interferometer and the AO system works for 
this limited inner region. 

For the two wavefronts, we assume the same unaberrated 
amplitude, 




and complex aberrations of similar level (variance) as 
(l A (x,y)) = (l E B (x,y)). (6) 

3.2. Unbalanced nulling interferometer (UNI) process 

The first effect of the UNI stage is a starlight attenuation by 
a nulling interference of the two beams. For the wavefront Eg 
(plane B) the modulus of the complex amplitude is slightly re- 
duced, for instance by an ND filter. At the output of this ND 
filter, the complex amplitude becomes, 

E w (x,y) = (l-g)(E°+s li (x,yj) (7) 

where g is a small fraction of unity, i.e. the unbalance factor. 
The beam is then phase-shifted, such that 

E B „(x,y) = j*(l-g)(E°+s B (x,yj) 

= -(l-g)(£° +eB (x,y)). (8) 

where we considered (f> — n. 

The two wavefronts E\(x,y) and E B "(x,y) are summed for 
destructive interference, so that at the UNI output plane (sub- 
scripted UNI), and using Eq.©, 

E um (x,y) = gE A +e A (x,y)-(l-g)s B (x,y), (9) 

where we omitted the reflection and transmission ratio of the 
beam combiner. 

There is one point we need to clarify here: the important ef- 
fect we want to achieve is to obtain an appropriate value for the 
modulus of the unaberrated component of the residual complex 
amplitude after the beam combination. In other words, it does 
not matter whether this residual is purely real or not. Therefore, 
in Eq. ©, the factor e I(S (l - g) can be complex, as long as we 
can control either g or <p (or both) to produce the desired un- 
aberrated modulus amplitude for the combined beam. In the 
following, we assume <p - n for simplicity. 
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We consider that defects of the optics for the unbalancing as 
well as for the phase shift optics are negligible, or included in 
£b ■ The aberrations of the beam combiner and the two optical 
paths in the UNI are treated in the same way using b a and s^- 

The unaberrated complex amplitude of Eq. (O and its in- 
tensity become 

= SE° A 

(10) 



0.4 



E° 

UNI 



(ID 



At this point, the starlight intensity is reduced by g 2 , e.g., if 
g = 0.1 the starlight is attenuated by 1/100. The complex aber- 
ration term is almost the sum of the two wavefronts as 

( euNi (x, y) = s A (x , y) - ( 1 - g) e B (x, y) 
U'uniM = (2 - 2g + /)(/£(*, y)), 

where Eq.© and the independence of e A (x,y) and es(x,y) are 
used for the latter equation. This situation corresponds to a 
"worst case" because if the aberrations are perfectly correlated 
so that e A (x,y) = s B (x,y), then I^ m (x,y) = g 2 I s A (x,y), mean- 
ing that the aberration term is well canceled as well as the un- 
aberrated amplitude. Therefore, we may actually get a better 
result in practice than the Eq. (fTTT i if the wavefronts partially 
correlate (see Appendix |B). We continue with the worst case 
formula proving the UN1-PAC process works well. 

In order to evaluate the dynamic range improvement by the 
UNI stage, the unaberrated and the aberrated component of the 
star should be compared with the planet. We introduce the in- 
tensity of the planet l p A , where we only consider the unaber- 
rated component, since its aberration component intensity is 
very small. We adopt the following expression for the planet 
intensity after the UNI as 



L \M\ 



a-ig+g 2 )! 1 : 



(12) 



by averaging for the constructive and the destructive conditions 
of the interferometer (see Appendix C). Then using Eqs. ( ITOl . 
(Q~TJl, and ( 112) . the dynamic range improvement through the 
UNI stage is expressed as 



7° 

J UNI 

L v,m 



1° 



v 

*UNI 



2-2g + g z 



(13) 



where the latter equation shows that the averaged aberration 
intensity is not affected by the UNI, while the star intensity gets 
closer to the planet intensity by a factor of g 2 /(2 - 2g + g 2 ). 

3.3. Wavefront aberration magnification by UNI stage 

The important phenomenon at the UNI stage we must empha- 
size is the wavefront aberration magnification. The normalized 
complex aberration after the UNI stage can be evaluated, using 
Eqs. dTOj and 0. 

guNifcy) _ 1 £A(x,y) - (1 -g) s B (x,y) 



^UNI 
^UNI 



—J 



£A(x,y) 



(14) 
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Fig. 4. Wavefront aberration magnification by the unbalanced 
nulling interference, (a) Complex amplitude of 400 points is 
simulated in a wavefront whose phase aberration rms is 0.0063 
radian and amplitude aberration rms is 0.004. (b) After the un- 
balanced nulling interference with a no-aberration wavefront 
with an amplitude ratio of 0.9, the phase aberration is easily 
seen as magnified, and (c) the amplitude aberration magnifica- 
tion is also obvious when the complex numbers are normalized 
by the unaberrated wavefront electric field, which is 0.1 E° in 
this case if we represent the original unaberrated field as E°. 
The magnified aberrations seen in panel (c) can be corrected 
again to the same quality as (a) by an AO system in the next 
stage. 

Here the UNI magnifies the normalized complex aberration by 
a factor 1 jg (recalling that g < 1) compared to the initial wave- 
front aberrations, or, in average by ^2-2g + g 2 /g. 

In order to obtain an appropriate magnification of the wave- 
front aberration, we chose a sufficiently large value of the un- 
balance factor g, so that EtjniOc, y) always has a positive inten- 
sity (modulus) without any phase singularities, |£'uni(-* : ,3 ; )I > 0. 
A little tighter constraint is also suitable, 



> |euNi(^,y)l 



(15) 



which means E\j^i(x, y) is distributed over half of the complex 
plane excluding the origin (see Pig . [3j ■ Within the condition of 
this equation, the smaller g, the more the star signal is attenu- 
ated and the aberration is magnified. 

An example of the aberration magnification by the UNI 
stage is shown in FigjU where one wavefront has /1/1000 
phase, 0.004 amplitude aberrations in rms, and an unbalance 
factory = 0.1. The initial aberration is \sA(x,y)/E A \ < 0.02 and 
the magnified one is \s\jNi(x,y)/E® Nl \ < 0.2 which will be com- 
pensated again to less than 0.02 by the PAC stage. Although 
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Fig. 5. Results of a simplified numerical simulation showing 
the improvement of the dynamic range in the image plane by 
the UNI-PAC method. The curves show the central star (bold 
line) and its speckle halo noise (solid line) normalized by the 
planet peak intensity (dashed line). Initially (left panel) the con- 
trast between the star and the planet is assumed to be 5xl0 9 
and the speckle halo lies 4x10 above the planet, which is a 
result of the phase aberration of /1/1000 rms and the amplitude 
aberration of 0.004 rms with a flat spectrum (30 x 30 square 
grids in the pupil). After the UNI is applied with the factor 
g — 0.1 (middle panel), the central star intensity is reduced by 
g 2 /(2 - 2g + g 2 )=0.0055 times. After the PAC process (right 
panel), the aberration-magnified wavefront is compensated to 
the same aberration rms as the input. The speckle halo noise is 
then decreased A:"=0.0055 times, that is, twice the planet peak 
intensity. The star profile is drawn as a point spread function 
of a non-modified pupil, which can be suppressed by a down- 
stream coronagraph. 



g = 0.02 can be adopted just to satisfy Eq.dTsb. a magnified 
wavefront aberration close to unity \eumix,y)/E9j m \ ~ 1 would 
be less favorable to perform the wavefront correction at the 
PAC stage. 

If g — > 0, the unaberrated amplitude tends to completely 
cancel out, E^ m — > 0, and the nulling result consists of only the 
aberration terms, £tjni(x, y) = sa(x, y) - eb(x, y) which are dis- 
tributed around the origin in the complex plane. In this case, the 
residual intensity in the pupil plane takes the form of a speck- 
led field with a small inhomogeneous amplitude and random 
phase fluctuations, resulting in the existence of zero amplitude 
points. Consequently, low light levels and zero intensity points 
with phase singularities make it difficult to measure and correct 
the wavefront with an AO system. This is exactly the problem 
situation we aim to avoid with the unbalanced nulling. 



3.4. Phase and amplitude correction (PAC) process 

After the wavefront aberrations have been magnified by the 
UNI stage, the goal is to compensate them to a level equivalent 
to SA(x,y)/E^ (the initial wavefront aberration) by a secondary 
AO system. In order to perform the phase and amplitude correc- 
tion (i.e., the PAC), the UNI is followed by a second AO system 
composed by two DMs (Fig.[TJ and a WFS which measures the 
wavefront shape and the intensity distribution. The compensa- 
tion level can be set differently according to the characteristics 
of each AO system. 



The PAC process is formalized by applying a wavefront- 
aberration-reduction function h(x,y) to the normalized com- 
plex aberration, 



E° 

^ PAC 

gpAcO,y) 



^UNI 



"PAC 



= k 2 



guNifey) 



E° 

^UNI 



(16) 



where A; is a kind of weighted average of the aberration- 
reduction function h(x, y) defined as 



k 2 



(\h(x,y) s VN1 (x,y)\ 2 } 
(leuNifey)! 2 ) 



(17) 



Since the residual aberrations before and after the UNI- 
PAC are caused by sensing, control, specifications of the DMs, 
or non-common path errors of each AO system, then spAc(x,y) 
or h(x,y) is not correlated with sum(x,y). 

After the PAC operation, the unaberrated complex ampli- 
tude can be 



E° 

11 PAC 

7° 

*PAC 



= I? 



(18) 



In the equations above, we assumed a negligible intensity 
change due to the amplitude control of the PAC DMs (under 
very small wavefront aberration conditions). If we were to con- 
sider this change in intensity, we should introduce a parameter 
^ such that Ep AC = i; E^ ]m with ^ ~ 1 in our equations. But its 
effect on the principle of the present concept should not be very 
large. In the following, we do not take this effect into consider- 
ation. 

In the PAC output plane, the complex aberration term be- 
comes 



(19) 



( s P Ac(x,y) = h(x,y) £ Vm (x,y) 
1('?acM = k 2 {l B VN1 (x,y)), 
In total the complex amplitude is 

EpAc(x,y) =gE° A + h(x,y) (sA(x,y) - (1 - g) e B (x,y)) . (20) 

For the planet, the intensity change caused by the PAC 
stage can be considered negligible for reasons similar to that 
for Eq.dH, 



f = f 



(21) 



Then the improvement of the dynamic range at the PAC stage 
is just the aberration component reduction expressed as, 



7° 

*PAC 
l PAC 



I P 



(22) 



L PAC 



L \M\ 



If we want to express the PAC stage performance compared to 
the initial AO-compensated wavefront level, we can write, 

2\ 



e P Ac(^,y) 



E° 

^PAC 



£a (x,y) 



E° 

A 



(23) 
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where r\ ~ 1 means that the residual wavefront aberration after 
the PAC is similar to that delivered by the first AO system. 
A smaller value of rj would indicate a better compensation by 
the PAC stage, and inversely. Then with Eqs. ( fT4b and ( fToT ). we 
obtain 



k 2 = rj 2 g 2 /(2-2g + g 2 ). 



(24) 



This expression shows that, when 77 ~ 1, A: 2 is a number less 
than unity almost completely defined by the unbalance factor g 
of the UNI stage. 

3.5. UNI-PAC total effect 

From equations (fT3l l and (l22l , the total reduction effect on the 
intensities of the unaberrated and aberrated wavefront compo- 
nents are respectively, 



7° 

1 PAC 
*PAC 



i p 



= k 



2-lg + g 1 



7 A 



(25) 



5 



Equivalently in the image plane, the central star profile attenu- 
ation factor is g 2 /(2-2g + g 2 ) and the average speckle halo at- 
tenuation factor is k 2 . These results can be seen in Fig. [5] based 
on a simplified numerical simulation where the corrected wave 
after the PAC stage is calculated by generating random residual 
aberrations. 

The very important conclusion is that the sequential com- 
bination of the UNI and PAC achieves an aberration compen- 
sation of the star wavefront amplitude beyond the correction 
capabilities of the individual AO systems used (including the 
non-common path errors). We wish to call it a virtual compen- 
sation. 

With this pre-optics, what a downstream coronagraph 
should do is to perform the rest of the total dynamic range, 
where the requirements for the wavefront aberration level and 
the star intensity suppression level are relaxed by the reduction 
factors of the unaberrated and the aberration wavefront (the star 
and the speckle) intensities in Eq. d25l l. The detailed contrast 
improvement in the image plane by this concept depends on 
the spectrum of the aberrations and the specifications of the 
AO system, although the (average) relative improvement of the 
dynamic range is the same as the pupil plane. 

4. Simulation of the UNI-PAC 

In this section we want to illustrate in a more global way 
how the UNI-PAC can perform, aside from purely analytical 
considerations. In our simulations, we use the Fourier optics 
approximation, and FFT to produce complex fields at conju- 
gated planes. The top row of Fig. [6] shows the distribution of 
the two complex amplitudes of the beams A and B. The ini- 
tial phase aberrations are uniformly distributed over pupil spa- 
tial frequencies (flat power spectrum) with an rms amplitude 
of A 1 1000, whereas the amplitude aberrations follow a /~ 25 
power-law with an amplitude of 1/200* of the average ampli- 
tude level. Beam B is shown after passing through a neutral 
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Fig. 6. Simulation of the UNI-PAC resulting in a ~ A/ 10000 
rms from input wavefronts quality of ~ ,4/ 1000 rms. The top 
row shows the complex amplitude distribution in the complex 
plane for the input Beam A, Beam B (including the unbalance 
factor g = 0.1), and for the unbalanced combined beam, i.e. 
at the UNI output plane. The middle row shows the amplitude 
(modulus) distribution of UNI output beam after the propaga- 
tion between DM1 and DM2, before (left), and after (right) the 
Gerchberg-Saxton algorithm convergence. DM1 compensates 
for amplitude fluctuations only, while DM2 compensates for 
the phase. The bottom row shows the PAC output wavefront 
complex distribution. Note that the real and imaginary axes 
scale has been reduced by a factor of 1/10. 



density filter with 8 1 % transmission, meaning that at the UNI 
output plane, the recombined beam intensity is reduced by a 
factor g 2 = 0.01, and the phase aberrations are amplified. The 
PAC stage uses two deformable mirrors (DM1 and DM2). DM1 
modulates the UNI output beam in phase only so that the cur- 
vature of the wavefront induces an amplitude modification by 
beam propagation to the location of DM2. Once the phase and 
amplitude after the UNI are measured, we use an FFT-based 
propagation algorithm combined with a recursive Gerchberg- 



Saxton (Gerchberg & Saxton 1972 and e.g. Fienup 1982 1 al- 
gorithm, in order to find the DM1 phase which converges to- 
ward a given pupil amplitude distribution (flat in this case). 
This step is shown in the middle row of Fig. [6] The left panel 
shows the initial propagated beam amplitude, and the right 
panel shows the same beam, but after convergence of the al- 
gorithm. DM2 compensates for the phase aberrations includ- 
ing those at the UNI output, as well as the effects of phase 
modulation by DM1. We did not use a realistic phase anal- 
ysis/correction scheme. We simply assumed that the sensing 
sensitivity was high enough to provide a relative phase quality 
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Fig. 7. Azimuthally averaged profiles of focal plane images at 
different locations in the UNI-PAC system: a reference Airy 
profile (dotted), the profile at the UNI (dashed) and the PAC 
(solid) planes. The thick dotted profile corresponds to the 
speckled halo (aberration component) before the UNI stage. 
The two thick lines at the bottom of this plot are (thick black) 
the image profile by the virtually-corrected PAC output wave- 
front with the ideal coronagraph, and (thick grey) the image 
computed by the initial wavefront with the artificial enhance- 
ment (i.e. with the complex aberrations reduced by a factor of 
1/10) combined with the ideal coronagraph. 

similar to that used for the input beams (Beam A and Beam B). 
In this simulation, the DMs are 32 x 32 actuators (793 within 
a circular pupil are effective over a total of 1024). The dis- 
tance between the two DMs in the PAC stage was set to a rela- 
tively arbitrary distance of 30 cm (propagation distance) while 
the collimated beam diameter was 1 cm. These parameters cer- 
tainly need to be better adjusted, especially regarding the DMs 
dynamical range, and the actuators influence functions, but this 
is not the purpose of this paper. The bottom of Fig.[6]shows the 
final corrected wavefront at the PAC output plane. It has ap- 
proximately the same wavefront aberration quality as the input 
beams, both for phase and amplitude. 

Figure[7] shows various plots corresponding to different 
steps of this UNI-PAC simulation. The top dotted line is the 
profile of the PSF computed from Beam A. For comparison 
purposes, the corresponding thick dotted curve is the profile of 
an image given by an "ideal" coronagraph (i.e., which perfectly 
cancels the unaberrated plane wave) with this input wavefront 
(around 10 7 contrast level). The dashed and continuous thin- 
ner curves are respectively the image profiles at the UNI out- 
put, and after the PAC correction. Note the small improve- 
ment of the PAC profile (continuous line) further away in the 
halo which means a better wavefront. Then, we applied the 
same "ideal" coronagraph to the PAC output, which is shown 
as the thick black continuous curve at the bottom of the fig- 
ure, where the virtual aberration correction improvement be- 
comes obvious. For comparison, we have generated a similar 
coronagraphic image, but from the Beam A data of the initial 
wavefront, with its aberrations artificially attenuated by a fac- 
tor of 1 / 10, i.e. a wavefront with a phase aberration of A/ 10000 



jj 9.97 - 

ft / 

9.96 - X 

o / 

* 9.95 I — 1 1 1 

450 500 550 600 650 

Wavelength (nm) 

Fig. 8. Wavelength dependence of the magnified aberration. 
The wavefront shape aberration after the UNI, ztjni> is shown 
for the input wavefront shape aberration za = 1 nm with zero 
amplitude aberration and g = 0.1. In the wavelength band of 
500 nm to 600 nm, the wavefront aberration is magnified about 
9.96 times with a small wavefront difference of about 0.01 nm, 
i.e., i/50000 level. 

and amplitude aberrations of 1/2000. This latter comparison 
demonstrates a potential of the present method, i.e. the virtual 
aberration correction improvement by the UNI-PAC is equiva- 
lent to the initial compensation. 

5. Discussion 

5.1. The virtual compensation 

In equation (1201 . at the PAC output plane, both the initial com- 
plex aberrations e\(x,y) and s%(x,y) are multiplied by a factor 
h(x,y). This can be interpreted as the result of an unbalanced 
nulling interference of two wavefronts written as 

f E A (x,y) = E°. + h(x,y) s A (x,y) 

(26) 

{E B ( X ,y) = El+h(x,y)eB( X ,y). 

In other words, it is equivalent to having a much better input 
wavefront (h times smaller), followed by an unbalanced nulling 
beam combination, but no further AO correction. At the UNI- 
PAC system output, one cannot distinguish whether the wave- 
front was very good at the input and inside the UNI optics or 
if it was processed by the UNI-PAC (see also Fig.[7]l. Thus one 
can consider that the present concept allows for a virtual com- 
pensation beyond the AO performance. 

5.2. Chromatism issues 

In this paper the formalism is constructed for one wavelength 
and the method is not perfectly achromatic. 

5.2.1. Magnification process 

If the wavefront shape expressed as z(x,y) (for instance, nm 
unit), is common to every wavelength because of mirror reflec- 
tions, then the phase aberration 6(x,y) = 2nz(x,y) /A is wave- 
length dependent. 

For the two beams, let us consider the wavefront shape 
aberrations of za (t-e., 8a = 2ttza M) and zb = (this latter 
value is set to zero for convenience, but it does not affect a 
more general result). Wavefronts before and after the UNI are 
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Fig. 9. Differential phase dependence of the magnified aber- 
ration. The wavefront shape aberration after the UNI, ztjni> is 
shown as a function of the additional phase change where the 
input wavefront shape aberration za — 1 nm with zero ampli- 
tude aberration under a condition of g — 0.1, B = 0.2 m, and 
A = 500 nm. Within the phase change of + 0.006 rad, the wave- 
front aberration is magnified about 9.95 times with a small 
wavefront difference of about +0. 1 nm, /1/5000 level. 

set to 



E A = E° + s 
E B = E° 
Evni = gE° + e, 



(27) 



so that the phase aberration magnification can be written (using 
the notation of AppendixlAl and the representation of Fig. lA.U 



tan 9 



UNI 



1 + e T /E° 

g + Sr/E 



tan 6U 



(28) 



where e r is the real part of e and #uni is the phase aberration 
after the UNI. The aberration magnification is made through 
the physical process in each wavelength and here we should 
check how the wavefront shape zuni are common in the wave- 
length band. The plot in Fig. [8] shows the wavefront aberration 
at a point in the pupil plane magnified by the UNI for an initial 
aberration za = 1 nm ( an d a zero amplitude aberration), for an 
unbalance factor g = 0.1. The differential wavefront shape af- 
ter the magnification is about 0.01 nm (/1/50000 level), within a 
wavelength band between 500 nm and 600 nm. Here the nonlin- 
earity of the tangent is the main cause of the wavelength depen- 
dence. The intensity-unbalanced nulling beam combiner has 
less dispersion, at least for an input aberration level of /1/1000 
with an aberration magnification factor of 10. 

5.2.2. Phase shifter and ND filter 

Another issue of the achromaticity is for the ji phase shifter and 
the ND filter in the UNI stage, i.e., a wavelength dependence 
of the factor -(1 - g) (e.g., (p and g) in Eq.®. 

Similar to the previous section, we consider a combination 
of two wavefronts by UNI written, 



E A = E° + s 

E B = e*E° 

E VN i = (l~(l~g)j*)E° 



(29) 



wherer t/c is a (wavelength-dependent) differential phase change 
in the wavefront B, then the aberration magnification is ex- 



pressed as 



£UNI 



E° 

UNI 



l-(l-g)e 



E 



(30) 



Here the magnification factor can be calculated for a given 
differential phase change independent of its origin (such as a 
wavelength dependence or an incident angle to an interferome- 
ter with a baseline described in the next section). The magnified 
wavefront aberration through the UNI process as a function of 
the differential phase change is calculated and shown in Fig|9] 
for the input wavefront aberration of 1 nm (/1/500) for phase 
and zero for amplitude. We find that the aberration is magnified 
to around 9.95 nm {A/50) with the enough small difference of 
about 0. 1 nm (/1/5000) within the differential phase change of 
±0.006 rad (±/l/1000). The wavelength dependence of the mag- 
nification described in the previous section is not more critical 
than this term. 

Note that for a given differential amplitude change writ- 
ten in a range of g = 0.1 + 0.001, i.e., +0.1% difference of the 
transmission of the ND filter in wavelengths, a similar level of 
the differential aberration magnification can be found. We leave 
a more detailed study of these differential magnification issues 
for a future work. 

5.3. Resolved stellar disk 

The formalism in Sect. [3] treated an on-axis point source as 
the central star. When the source has a resolved diameter, 
we should consider off-axis rays where the wavefront aber- 
rations are magnified by different factors by the UNI stage. 
Considering the baseline B between the two interferometer 
apertures, the two wavefronts have an additional phase differ- 
ence of iff — 2tt& smx I A, where x is the incident angle of the 
off-axis light. So that we can refer the previous subsection to 
address this point. 

For instance, with a baseline of 0.2 m (e.g., 
Sha o et al. 20041 > 0.5 mas tilted wavefront which comes 
from the limb position of the Sun at lOpc has a phase dif- 
ference of t/f = 6.1xl0~ 3 radian (/1/1000) between the two 
apertures for the wavelength of 500 nm. Therefore the range of 
the incident anlgle for the Sun at lOpc almost correeponds to 
the phase range shown by the curve in Fig |9] with the enough 
small differential magnification of about /1/5000. This means 
that the stellar diameter might not affect the performance 
for a short baseline interferometer under small aberration 
conditions. 



5.4. Larger phase aberrations 

Even for a large wavefront phase aberration of A/20 and am- 
plitude aberration of 0.1, the present method is still effective 
with the factor g = 0.3, where the central star extinction and 
speckle intensity reduction becomes about 0.05. The aberration 
level is similar to a ground-based AO condition and the gain of 
0.05 is essential for exo-planet observations from the ground, 
e.g., extending the observable range from 10~ 5 to 5xl0~ 7 . If 
the PAC AO reaches higher bandwidth than the telescope-side 
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AO system (e.g., |Guyon 2006) 1, it might be possible to apply 
the UNI-PAC method to ground-based observations. 

6. Conclusions 

We have presented the principle of a pre-optics concept for pre- 
cise wavefront aberration reduction in front of a coronagraph 
in a terrestrial-planet finding telescope. It consists of a com- 
bination of an unbalanced nulling interferometer and a two- 
deformable mirror AO system with a conventional pupil plane 
wavefront sensor, where it does not meet the low intensity prob- 
lem for wavefront measurements. This method reduces both the 
source intensity as a nulling coronagraph and the (speckle halo) 
noise intensity produced by wavefront aberrations. 

The wavefront magnification phenomenon at the unbal- 
anced nulling interference makes it possible to correct the 
wavefront precisely beyond the capabilities of employed AO 
systems with non-common path errors. In reality, by using 
the UNI-PAC, the specification for all of the coronagraph op- 
tics can be relaxed which is probably key to developing a 
cost-effective exo-planet detection system, although the per- 
formance of the whole system should be investigated for each 
case. 

A candidate instrument where the UNI-PAC system can be 
used would be a modified Michelson interferometer in a col- 
limated beam of an off-axis telescope (IShao et al. 2004t . We 
also expect the present concept to be applicable with several 
other coronagraphic concepts instead of the UNI stage by a 
Michelson beam combiner, i.e., with more general contrast- 
reduction setups (see for example lAbe et al. 20061 1. but we 
leave this discussion for a future study. 

Some other considerations, such as achromaticity, resolved 
star nulling, planet image quality, and AO performance limita- 
tions should be investigated. 

Acknowledgements. LA is supported by Grant-in-Aid (Nos. 
160772048007 and 160871018002) from the Ministry of Education, 
Culture, Sports, Science, and Technology (MEXT) of Japan. JN 
and NM are supported by Grant-in-Aid (No. 19360036 and the 
Priority Areas "Development of Extra-solar Planetary Science") from 
the MEXT. The authors are grateful to Profs. M. Yoshizawa and 
M. Tamura for their constant encouragements with this work. 

Appendix A: Wavefront aberration expressions 

The complex amplitude of an electric field in a pupil plane with 
a "classical" wavefront phase aberration 9{x,y) and an ampli- 
tude aberration a(x, y) can be written as 



E(x,y) = E°(l + a(x,y)) e Wt -- v) 



(A.l) 



Here E° is the unaberrated amplitude (usually assumed as a 
real number), and the effect of the aberrations can be defined 
as a deviation e between E and E° noted, 



s{x, y) = E(x, y) - E , 



(A.2) 




Fig. A.l. Correspondence of the electric field definitions. 

e and the "classical" wavefront aberrations can be seen in 
Fig lA.ll written as 



tan 6(x,y) = - S ^ E ljAI 

1 + £r/£° (A.3) 

a(x,y) = V(l +s r /E ) 2 + ( Ei /E ) 2 - 1 , 
and 

sJE° = (l + a(x,yfj sin 9(x,y) 

s t /E° = (l +a(x,y)) cos0(x,y) - 1 (A.4) 

s(x, y) = s r + i £j 

where i is the imaginary unit. e r and e; are real numbers which 
indicate the real and the imaginary part of s{x, y), respectively, 
and where (x, y) is not indicated for simplicity. 

E° can be defined as a spatial ensemble average of E within 
the pupil plane, E° = (E(x, y)) so that the aberration component 
has a zero mean, (e(x,y)) = 0, whereas (a(x,y)) and (6(x,y)) 
are not always necessarily zero. Under very small wavefront 
aberration conditions \s/E a \ <k 1, as for the present method, 
approximations 



6{x,y) ~ Ei/E° 
a(x,y) ~ Sy/E 



(A.5) 



can be used. Therefore the normalized complex aberration, 
s(x,y)/E° is a good indicator of the wavefront aberrations for 
our concept. 

Appendix B: Correlation between the two 
wavefront aberrations at the UNI stage 

At the UNI stage, if the complex aberrations of the two wave- 
fronts strongly correlate and have almost the same amplitude, 
£A(x,y) ~ es{x,y), the difference between the two wavefronts 
at planes A and B tends to zero. Then the aberration component 
of the UNI output becomes 



(B.l) 



which is the complex aberration component in Eq.([TJ, at the 
beginning of the formalism in Sect[3] Relationships between 



( e vm (x,y) ~ gs A (x,y) 

IfeiM ~ 8 2 (l s A (x,y)), 

instead of Eq. (fTTT i. This means that a simple nulling inter- 
ference can reduce both the star intensity and the speckled 
halo. This constitutes a better situation for achieving a deeper 
nulling without taking being concerned with the aberration, but 
it would be unrealistic. 

If both wavefront aberrations are perfectly uncorrelated, 
the unbalanced nulling combination produces a larger complex 
aberration term. This is the "worst case" situation expressed by 
Eq.Cjj)- 
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When the two wavefronts are taken at different positions in 
a telescope primary or a telescope-side AO system, they may 
partially correlate and the nulling combination will provide re- 
sults intermediate between those described above. 

Appendix C: Planet intensity after the UNI stage 

We examine how an off-axis object (a planet) is affected by 
the UNI stage. We are working in an interferometric beam 
combination scheme (two sub-pupils extracted from a mono- 
lithic telescope mirror, or using diluted apertures). During a 
"blind search" observation, we must consider the rotation of 
the nulling interferometer to explore the planet in every direc- 
tion, which results in an averaging of the planet signal between 
the (pseudo-)nulled and constructive conditions of the interfer- 
ometer for a given baseline B (e.g., Bracewell 1978). Here the 
planet should also be sufficiently off-axis, i.e. more than ~ A/B 
to prevent the planet staying in the central nulled region. 

We introduce the complex amplitude of the planet 
E' A (=£g) and its intensity l p A , where we only consider the un- 
aberrated amplitude, and neglect the aberration term, since its 
intensity is very small. The average of the planet intensities 
for the constructive, \E A - (1 -g)E P ] \ 2 , and the destructive, 
| E A + (1 — g) £g | 2 , conditions after the UNI becomes 

I P m = (2-2g + g 2 )I P A , (C.l) 

which should be used for calculation of the dynamic range im- 
provement. 
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